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Summary 

A simple modification of Gilbert’s (1970) formula to account for slight 
lateral heterogeneity of the Earth leads to a convenient formula to calculate 
synthetic long-period seismograms. Partial derivatives are easily calculated, 
thus the formula is suitable for direct inversion of seismograms for lateral 
heterogeneity of the Earth. 

Derivation 

Gilbert(l970), using the classic results due to Rayleigh and Routh, 
derived a convenient formula to calculate synthetic seismograms by normal 
mode summation. We show that a simple modification of his approach to account 
for slight lateral heterogeneity is possible. 

The equation of motion is given by 

(p 0 +6p)d?u = (H„+H)u+f(t,x,), (1) 

where p 0 is density and H a is an operator appropriate in a spherically sym- 
metric Earth, while 6p and H are their pejrturbations. The source is described 
by / (t ,x s ). The Laplace transformation with respect to time yields 


(p 0 +6p)p z u = (Ho +H)u i\f (p ,x s ), 


( 2 ) 
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where we used the bar to denote the transformed quantities. 

Now we use the eigenfunctions of the spherically symmetric Earth, iZ^ 0 ) 
(n=l,2,...), to expand u, where satisfies 

~PoU%ui 0) = H 0 u£°\ (3) 

The normalization of is done by 

f p.f&'ajPdV = «5 mn . (4) 

E 

where * means complex conjugation. We express u by 

G = (5) 

n 

Substitute (5) in (2), multiply and integrate over the whole volume of 
the Earth. We obtain 

(p 2 +fc>n)on = S a s(< n [ ^ | s >— J) S <7X 1 <5p I s >)+/ 7 W . (6) 


where 


<n 1 6p | s > = f dpujjfi u^dV, 

E 

<n I H ] s > = fu^^Hu^dV, (7) 

E 

K = fun 0) *f(p,x s )dV. 

E 

Considering that the first two terms on the right-hand side of (6) are small, the 
first approximation of a ^ is given by 
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Then the next higher approximation is obviously 


a** = E(<*l*l.>V<.l«Pl.»5 5C ^s 5gr + 


O) 


We use in this order of approximation in the following derivation. 


Formula for the step function type source 

If wo assume the step function type source, F n should be replaced by F n /p> 


a* = £(<n|/i | s >—p 2 <n | <Sp | s >) 


K 


K 


p(p 2 +Us)(p 2 +U%) p(p Z +U%)' 


( 10 ) 


For a point source described by a second order tensor JSfy , F n is expressed by 


F'n ~ XI ^ V * 

a 


where ey is the strain of a mode at the source, as is well known, 

Substituting (10) in (5) and inverting to the time domain, we obtain 
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<n | // j s > - --—sin w n t ( <n | H |s>+*c^<n | 6p\s>) 
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cos<o*f (<n| jy r |s>+w^<n |<5pjs>)) 
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Note that the first term is equivalent to Gilbert’s (1970) formula and the rest are 
corrections due to the heterogeneity of the Earth. 

Simplification 

A few simplifications to the formula (11) are possible. First of all, static 
terms are not necessary in many applications, since they are usually filtered out 
in the real data. Secondly, if we assume that u n = o a occurs only within the 
same multiplet, then the sum over s is replaced by the sum over the azimuthal 
order number m only. Also for a short time approximation, the secular term can 
be absorbed in cosw n £ as 


cosw n f = cosw n f -f — — sinci> n t 


v(\H\m.>+u%< iu\6p\m>)F„ 
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where 


EE“K>(< m\H\m>+ul< m\6p\m>)F m 

= (13) 

E aSf. 
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and it is understood that the sum over m and rri are taken within the multiplet 
specified by n. Note that the last term in (13) devided by 2w n , is the multiplet 
location parameter (Jordan, 1970, Woodhouse and Girnius, 1981). Thirdly, the 
last term in (ll) can be exchanged with the last term in the case of u ^ for com- 
putational convenience. Then we have 
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to the accuracy of the above assumptions, s and n specify multiplets and m 

and mf refer to certain singlets within each multiple! respectively. The 

definitions of the new symbols should be clear in the context. 

It is straight-forward to use (14), since some relevant formulae for <n j/7js> 
and <n]<5p|s> are available in the literature (e.g. Woodhouse and Dahlen, 1978) 

and others are not so difficult to derive. Also it is possible to express <n\H\s> 

in terms of elastic constants, <5\ and 6ja, thus (14) is useful to invert the data for 
6p, 6X and 6p,. 

Care must be taken in (14), because (14) breaks down if co s is very close to u n . 
This is because (14) assumes that the terms from summation over s are small 
perturbations. Such a case can be treated in (ll), however, by incorporating 
those modes in the case for v n — v s . 

Lastly, the equation (14) seems to be equivalent to the formula in Wood- 
house(l983), although the two derivations are very different. 
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